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Abstract 

In 1947 Mills proved that there exists a constant A such that Lj4 3 "J is a prime 
for every positive integer n. Determining A requires determining an effective Hoheisel 
type result on the primes in short intervals — though most books ignore this difficulty. 
Under the Riemann Hypothesis, we show that there exists at least one prime between 
every pair of consecutive cubes and determine (given RH) that the least possible value 
of Mills' constant A does begin with 1.3063778838. We calculate this value to 6850 
decimal places by determining the associated primes to over 6000 digits and probable 
primes (PRPs) to over 60000 digits. We also apply the Cramer-Granville Conjecture 
to Honaker's problem in a related context. 

1 Introduction 

In 1947 Mills [19] proved that there exists a constant A such that L^ 3 "J is a prime for every 
positive integer n. Mills' proof used Ingham's result [17] that there is always a prime in 
the interval (x, x + kx 5 ^ 8 ) for some constant k. Since k was not determined explicitly, Mills 
could not determine A or even a range for A. Many authors citing Mills' constant either 
follow his example and remain silent about A's value (e.g., [11], [24], [31]), or explicitly state 
that it is unknown (e.g., [1], [5]). But recently a few authors have begun to state that A 
is approximately 1.3063778838 (e.g., [8], [10], [28], [33] and Sloane A051254). They do this 



despite the fact that it is currently impossible to compute A without any further unproven 
assumptions, and despite the fact that for each c > 2.106, there are uncountably many 
possible values A for which [A c ™ \ is prime for each positive integer n. 

The authors who approximate Mills' constant A also implicitly define Mills ' constant to 
be the least A such that L^ 3 ' ] is prime for all positive integers n. Mills' original article 
contained no numerics — it only showed that such an A exists. 

In Section 2 we observe that, assuming the Riemann Hypothesis, it easily follows that 
there is at least one prime between x 3 and (x + l) 3 for every x > 0.26 (Lemma 5). We then 
use this to determine Mills' constant A to over 6850 decimal places in Section 3. This requires 
a substantial computing effort to calculate the associated Mills' primes [A 3 ™ \ explicitly for 
n = 1, 2 . . . , 10. In particular we find prove the following: 

Theorem 1. Assume the Riemann Hypothesis. The minimum Mills 7 constant (for the ex- 
ponent c=3) begins with the following 600 digits: 



1.3063778838 
3644368053 
0729696093 
1943534578 
2361522359 
5279707992 
3310565215 
0405320517 
1227499787 
6091500412 
0736114409 
6154666880 



6308069046 
7599664340 
8103086882 
7110033188 
0622186016 
5631721527 
3513186684 
6890260257 
6664385157 
3788524793 
4137650364 
4874760951 



8614492602 
5376682659 
2388614478 
1405093575 
1085667905 
8412371307 
1550790793 
9344300869 
6619143877 
6260880466 
3793012676 
4410790754 



6057129167 
8821501403 
1635348688 
3558319326 
7215197976 
6584911245 
7238592335 
5290636205 
2844982077 
8815406437 
7211713103 
0698417260 



8458515671 
7011973957 
7133922146 
4801721383 
0951619929 
6317518426 
2208421842 
6989687262 
5905648255 
4425340131 
0265228386 
3473107746 



We also find the two probable-primes which are (most likely) the next terms in the sequence 
of Mills' primes. The last of these has 61684 digits. 

Concerning other problems that involve prime gaps, let po = 2, p n be the n-th odd prime 
number and g n = p n +i — Pn, the gap between consecutive prime numbers. The Cramer- 
Granville conjecture [16] states that g n < M\og 2 n for some constant M > 1. The last 
section is an application of this conjecture on prime gaps to Honaker's problem. 

Honaker's problem [9] is to find all trios of consecutive prime numbers p < q < r, 
such that p\(qr + 1). We call such triples of consecutive prime numbers Honaker trios and 
conjecture that there are only three Honaker trios: (2,3,5), (3,5,7), and (61,67,71). We 
establish that there are only three Honaker trios (p, q, r) for p < 2 x 10 17 and prove: 

Theorem 2. The Cramer- Granville conjecture with any constant M implies that there are 
only a finite number of Honaker trios. If also M < 199262, then there are exactly three trios. 
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2 Gaps Between Prime Numbers 

Let s(q) represent the next prime after q: s(q) = min{p prime | p > q}. Define the sequence 
of maximal gap primes as follows: 

( 2, if ./ = 0; 

Qj = < 3, if j = 1; 

[ min{p prime | s(p) — p > s(qj- 1 ) — qj-i}, if j > 1. 

In other words {qk} is the sequence of primes that are followed by maximal gaps (Sloan 
A002386). This sequences starts 2, 3, 7, 23, 89, 113, 523, 887, 1129, 1327 ... and exists for 
all subscripts k because limsup p _> +00 s(p) —p = +oo [15]. The list of maximal gap primes 
qj has been extended through all primes below 2 x 10 17 by Herzog and Silva [30] with much 
of their work verified by Nicely and others [22, 23]. 
With these definitions it is clear that 

sjp) ~ V < s(q k ) - q k 
\og 2 p ~ \og 2 q k 

for every q k < p < q k +\- Now using Nicely's tables of maximal gaps [22] we can easily verify 
the following result: 

Lemma 3. For 11 < p n < 2 x 10 17 , we have < 0.92064. 

— fn i log p — 

To explicitly determine Mills' constant we must know that there is always a prime between 
each pair of successive cubes from some established point onward. This follows easily from 
a result of Schoenfeld [29]. Recall 

li(x) := lim I h 



Lemma 4 (Schoenfeld). Assume the Riemann Hypothesis. For x > 2657, we have 

y/x\0gX . y/x\0gX 

H x ) a — < n x ) < H x ) + — a — • 

8/T Sir 

(Much stronger bounds are now available [27], but this is sufficient for our current needs.) 
Using this Lemma we see that if x > 2657 1//3 , then 

ir((x + l) 3 ) - 7r(a; 3 ) > li((x + l) 3 ) - li(x 3 ) - ^-{x + 1) 3/2 log(x + 1) 

47T 

> / ^---^-(x + l) 3 / 2 log(x+l) 
J x z \ogt 4vr v 

> —— — — (x + 1 plog x+l . 

31og(x) 47T 

This last lower bound is an increasing function of x and it is greater than one for x > 2657 1 / 3 . 
After using a computer program to check the smaller values of x we have shown the following: 
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Lemma 5. Assume the Riemann Hypothesis. There is at least one prime between x 3 and 
(x + l) 3 for every x > Xq = 2 1 / 3 — 1. 

This lemma states that, under the Riemann Hypothesis, there are prime numbers between 
consecutive cubes; we will see that this is necessary to calculate Mills' constant in the next 
section. 

Without the Riemann Hypothesis it is still possible to use explicit upper bounds on the 
zeta function [6, 13] to get a version of Lemma 5 with a far larger value of xq, roughly 
10 60ooooooooooooooooo 26]. In the next section we will show that calculating Mills' constant 
involves explicitly finding a prime larger than the cube-root of this bound. This bound 
so dramatically exceeds current computing abilities, that any current calculation of Mills' 
constant must involve an unproven assumption. This bound xq should get much smaller as 
the bounds on the zeta function are improved. 



3 Mills' constant 

We begin this section with a lemma. 

Lemma 6. If x > 1 and c > 2, then 1 + x c + x c ~ l < (1 + x) c . 

Proof. Dividing by x c and replacing x with 1/x we arrive at the equivalent inequality < 
(1 + x) c — (1 + x + x c ) (0 < x < 1). The inequality clearly holds when c = 2 (because it 
reduces to x > 0) and when x — 0. Now if x > 0, differentiate the right side with respect 
to c to get (1 + x) c log(l + x) — x c \og{x), which is clearly positive, so the inequality above 
holds for all c> 2. □ 

It will be useful to next recall a proof of a simple generalization of Mills' theorem. 

Theorem 7. Let S = {a n } be any sequence of integers satisfying the following property: there 
exist real numbers Xn and w with < w < 1, for which the open interval (x, x + x w ) contains 

there is a number A for which \_A C J is a subsequence of S. 

Proof. (We follow Ellison & Ellison [12].) Define a subsequence b n of S recursively by 

(a) b\ is equal to the least member of S for which b\ is greater than xq. 

(b) b n+ i is the least member of S satisfying b c n < b n+ i < b c n + b™°. 
Because c > 1/(1 — w) and c > 2, (b) can be written as: 

b c n < b n+l < i + b n+1 < i + b c n + bz < i + b c n + b c ~ l < (i + b n y 

the last inequality following from Lemma 6. For all positive integers n we can raise this to 
the c~( n+1 Hh power to get 

b c n < b c n+l < (i + b n+l y < (i + b n y . 
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This shows that the sequence-^ } is monotonic and bounded, therefore converges. Call 
its limit A. Finally, 

b n < A° n < bn + 1, 

so [A cn J = b n , yielding the chosen subsequence of S and completing the proof. □ 

When constructing the sequence {b n } in the previous proof, the condition (a) can be 
relaxed in two important ways. First, by replacing the words 'the least' with 'any', we see 
there are infinitely many choices for b±, and hence for the resulting value A. Second, it is not 
necessary that b\ > xq be satisfied by b±, as long as the terms satisfying (b) exist to a term 
b n which does satisfy b c n > xq. This will someday be important in removing the assumption 
of the Riemann Hypothesis from Theorem 1 by calculating a sequence of primes {&«}™ =1 that 
extends to the cube root of the bound x® discussed at the end of the previous section. 

It is proved by Baker et. al. [3] that p n +i ~Pn = 0(p^ 525 )- From the above theorem, one 
obtains the following proposition. 

Proposition 8. For every c > 2.106, there exist infinitely many A's such that [A ™ \ is a 
prime for every n. 

Wright [32] showed that the set of possible values of the constants A and c in this 
proposition (and several generalizations [25], [28], [31]) have the same cardinality as the 
continuum and are nowhere dense. So authors approximating "Mills' constant" must first 
decide how to choose just one. Mills specified the exponent c = 3. Mills also used the lower 
bound k 8 for the first prime in the sequence, where k is the integer constant in Ingham's 
result [17]. So to follow his proof literally would require that the first prime be 257 (and 
producing a constant A m 6.357861928837). But all authors offering an approximation agree 
implicitly on starting with the prime 2 and then choosing the least possible prime at each 
step, so we take this as the definition of Mills' constant. 

The sequence of minimal primes satisfying the criteria of the proof with c = 3 (Sloan 
A051254) begins with 

bi =2, 

b 2 =11, 

b 3 =1361, 

h =2521008887, 

b 5 =16022236204009818131831320183, 
b 6 =411310114921510480003052953791595317048613\ 
9623539759933135949994882770404074832568499. 

These first six terms were well known. We have now shown that assuming the Riemann 
Hypothesis, prime numbers exist between consecutive cubes, so (with RH) we know this 
sequence can be continued indefinitely. 

To make these fast growing primes b n easier to present, define a sequence a n by b n+ i = 
bl + a n . The sequence a n begins with 3, 30, 6, 80, 12, 450, 894, 3636, 70756 (Sloan A108739). 
The primality of the new terms bj, b$ an d bg (2285 digits) were proven using the program 
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Titanix [18] which is based on an elliptic curve test of Atkin ([2], [20]). The test for 6 g , which 
at the time was the third largest proven 'general' prime, was completed by Bouk de Water 
in 2000 using approximately five weeks of CPU time. The certificate was then verified using 
Jim Fougeron's program Cert-Val. 

Frangois Morain has verified the primality of the next term b\o (6854 digits, July 2005) 
using his current implementation of fastECPP [14, 21]. The computation was done on a 
cluster of six Xeon biprocessors at 2.6 GHz. Cumulated CPU time was approximately 68 
days (56 for the DOWNRUN, 12 for the proving part). 

In late 2004 Phil Carmody used his self-optimizing sieve generator to generate an appro- 
priate sieve. Using it and pfgw he verified the above results (as PRPs) and found the next 
two probable-primes in the sequence b n . They are the values corresponding to a w = 97220 
and an = 66768. These yield probable primes of 20562 and 61684 digits respectively, so 
their primality may not be proven for some time. 

The first ten terms of {b n } are sufficient to determine Mills' constant A to over 6850 
decimal places because we have from the proof of Theorem 7: 

bC <A<(b n + iy- n . 

A quick calculation now completes the proof of Theorem 1. 



4 Honaker's Problem 

Let g{p) be the length of the prime gap after p: g{p) — s{p) — p. 

Cramer's original conjecture says that g(p) = (1 + o(l))log 2 p. This conjecture could 
be too strong to be valid, so Granville [16] proposed there is a bound M for which g{p) < 
M log 2 p (we refer to this conjecture as the Cramer- Granville conjecture). As we discussed 
in Section 2, M = 0.92064 is valid for 11 < p < 2 x 10 17 . However, the constant involved in 
this conjecture may be M = 2e~ 7 ~ 1.123 [16]. 

Proof of Theorem 2. Searching for Honaker trios with p < 50 we find only (2, 3, 5) and 
(3, 5, 7). In what follows we will assume (p, q, r) is a Honaker trio and p > 50. 

Reformulate the problem by letting q = p + 2k and r = q + 21. Since p\{qr + 1), we have 
p\(Ak 2 + Akl + l). Thus, 

p<Ak 2 + Akl + l. (1) 

If k > I, then p < 8k 2 + 1 and g(p) = 2k > ^{p - l)/2. Otherwise k < I, p < 
4(Z - l) 2 + 4(Z - 1)/ + 1 and again g(q) = 21 > y/(p- l)/2. These gaps are bounded by the 
Cramer- Granville conjecture, so either p or q must satisfy: 

M > \/^V^ (2) 
V 2 log p 

The function on the right is increasing for p > 50, and is unbounded, so for any fixed M the 
number of Honaker trios is finite. 
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Next suppose p < 2 x 10 17 . By Lemma 3 we know M = 0.92064 will suffice in this 
range, so equation (1) gives p < 14142. A computer search in this range adds the third trio 
(61, 67, 71). Thus there are exactly three trios with p < 2 x 10 17 . 

Finally, since any additional solution must satisfy p > 2 x 10 17 , equation 2 shows M > 
199262. This completes the proof of Theorem 3. □ 

Clearly this simple proof can be extended to make similar statements about finding k 
consecutive primes for which one of the k divides the product of the other k — 1 plus or 
minus a fixed integer. 
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